Mon. Not. R. Astron. Soc. 000, [TP?? (2012) Printed 3 March 2013 (MN IATeX style file v2.2) 



On the Cooling Trend of SGR 0526-66 



Tolga Giiver 1 , Ersin Gogiis, 1 and Feryal Ozel 2 

1 Sabanci University, Faculty of Engineering and Natural Sciences, Orhanh Tuzla 34956 Istanbul Turkey 
£N| 2 Department of Astronomy, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721 

o 

(N 

, ^ ' ABSTRACT 

We present a systematic analysis of all archival Chandra observations of the 
soft-gamma repeater SGR 0526—66. Our results show that the X-ray flux of 
SGR 0526-66 decayed by about 20% between 2000 and 2009. We employ physi- 
cally motivated X-ray spectral models and determine the effective temperature and 
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the strength of the magnetic field at the surface as kT = 0.354 



+0.031 
-0.024 



keV and 



B = (3.73^Q'Jg) x 10 14 G, respectively. We find that the effective temperature remains 
constant within the statistical uncertainties and attribute the decrease in the source 
flux to a decrease in the emitting radius. We also perform timing analysis to measure 
the evolution of the spin period and the period derivative over the nine year interval. 
We find a period derivative of P = (4.0 ±0.5) xl0~ u s s , which allows us to infer 
the dipole magnetic field strength and compare it with the one determined spectro- 
scopically. Finally, we compare the effective temperature of SGR 0526—66 with the 
expected cooling trends from magnetized neutron stars and suggest an initial magnetic 
field strength of 10 15-16 G for the source. 
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1 INTRODUCTION 

Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pul- 
sars (AXPs) are the prime representatives of magnetars - 
a class of neutron stars that are thought to be powered 
by the decay of their superstrong magnetic fields (B~10 14 - 
10 15 G, Duncan & Thompson 1992). Aside from emitting 
energetic bursts in hard X-rays and soft gamma rays, magne- 
tars are bright X-ray sources (L^IO^-IO 36 erg s _1 ), emit- 
ting pulsed X-rays either persistently or episodically (see 
Woods & Thompson 2006 and Mereghetti 2008 for detailed 
reviews). X-ray observations of most of SGRs and AXPs over 
the last thirty years have shown that these sources are not 
steady emitters, exhibiting flux variability generally in con- 
junction with bursting activities (see e.g., Rea & Esposito 
2011). 

All known SGRs have been discovered when they went 
into a bursting phase, during which they emit repeated en- 
ergetic bursts in soft gamma-rays. SGR 0526—66 was not 
an exception: it was discovered as the first SGR with an 
extremely bright gamma-ray burst on 1979 March 5 that 
was followed by a tail emission clearly pulsating at 8.1 s 
(Mazets et al. 1979). The burst location is within the super- 
nova remnant, N49 in the Large Magellanic Cloud (Evans 
et al. 1980). The source exhibited 16 more bursts until 1983 
April 5 with much lower intensity than that of the March 
5th event (Aptekar et al. 2002). Since then, no bursts have 
been detected from SGR 0526-66. 



Due to its large distance and the fact that the source is 
embedded in a supernova remnant, spectral and timing stud- 
ies of SGR 0526—66 have been limited. The persistent X-ray 
source was identified with ROSAT, yielding the earliest view 
of the source as it is in burst quiescence (Rothschild, Kulka- 
rni & Lingenfelter 1994). SGR 0526-66 has been the target 
of two Chandra X-ray Observatory (CXO) observations in 
2000 and 2001. Thanks to the superb angular resolution of 
CXO, X-ray spectroscopy and timing investigations of the 
persistent source became possible. Based on these observa- 
tions, spin period of « 8.04 s was measured and a tentative 
P of « 6.5x10"" s s" 1 was inferred (Kulkarni et al. 2003). 
X-ray spectrum of the SGR 0526—66 could be described well 
with the empirical model of an absorbed blackbody plus a 
power law component. The resulting power law index was 
much steeper than that of other SGRs, while it was similar 
to the power law indices of AXPs, which led Kulkarni et al. 
(2003) to suggest that SGR 0526—66 is in a transition from 
the SGR-like phase to an AXP-like phase. 

SGR 0526—66 and its associated supernova remnant, 
N49 was recently observed with XMM-Newton. Tiengo et al. 
(2009) reported that the X-ray spectrum of the source does 
not show any significant changes compared to XMM-Newton 
observations performed in 2000 and 2001 and SGR 0526—66 
is still a bright persistent source, emitting at an X-ray lu- 
minosity of w 4xl0 35 erg s . Very recently, Park et al. 
(2012) provided an in depth analysis of the supernova rem- 
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Table 1. Chandra Observations of the SGR 0526-66. 



Date 


Observation 


Exposure 




ID 


(ks) 


2000-01-04 


747 


43.9 


2001-08-31 


1957 


53.4 


2009-07-18 


10123 


28.2 


2009-07-31 


10808 


30.2 


2009-09-16 


10807 


27.3 


2009-09-19 


10806 


27.9 



nant using four Chandra observations that were obtained 
in 2009. They refined the Sedov age of N49 as r Sc d « 4800 
yr. Park et al. (2012) also analysed the X-ray spectra of 
SGR 0526—66 using phenomenological models, such as the 
sum of a blackbody and a power-law or two blackbody mod- 
els. They find that the X-ray flux of SGR 0526-66 var- 
ied by about 15 % between the observations performed in 
2000/2001 and 2009. Further including archival ROSAT ob- 
servations they detect a decay in X-ray flux by about « 20- 
30% for the last 17 years. 

Here, we report on our systematic analysis of all archival 
Chandra observations of SGR 0526—66 to unveil the X- 
ray spectral characteristics of the source when it is in deep 
burst quiescence. We model the X-ray spectra with the 
Surface Thermal Emission and Magnetospheric Scattering 
(STEMS) model (Giiver et al. 2007, 2008, 2011), to deter- 
mine the effective surface temperature, the surface magnetic 
field strength, and the long term X-ray flux of SGR 0526—66. 
We also perform timing analysis to construct the spin period 
evolution of the source and obtain the inferred dipole mag- 
netic field strength. Finally, we compare our results with the 
expected cooling trends from magnetized neutron stars. 



2 OBSERVATIONS AND DATA ANALYSIS 

We list the observations used in our study in Table [T] Note 
that there are two more Chandra observations that contain 
SGR 0526—66 in the field of view. However, the source was 
observed off-axis by about 4" in one of the observations (Ob- 
sID: 1041). The effects of vignetting do not allow a reliable 
source extraction without any contribution from the super- 
nova remnant. The other observation (ObsID: 2515) has only 
7 ks of effective exposure, that is too short to provide high 
enough signal-to-noise X-ray spectrum. Therefore, we ex- 
clude these two pointings from our investigations. 

We reprocessed all of the observations with the CIAO 
software suit, version 4.2 and CALDB 4.3.0 using the re- 
cently introduced chandrajrepr^ tool. This tool automates 
the creation of new bad pixel file and level 2 event file. We 
applied barycentric correction to each event file using the 
axbary tool. To extract source spectra and lightcurves we 
selected a circular region centering the neutron star with a 
radius of 2" , as shown in Figure [T] For the background we 
used an annular region that is centered on the coordinates 
of the neutron star and covers the range starting from 2.5 
to 5"from the center. Finally, we selected data from a larger 

1 http:/ /cxc. harvard.edu/ciao/ahelp/chandra_repro. html 



annular region: centered on the neutron star and covering 
the range starting from 4 to 11", to generate an X-ray spec- 
trum of the SNR and obtain an independent measurement of 
the hydrogen column density, as detailed below. Source and 
SNR extraction regions are shown in Figure [1] We extracted 
X-ray spectra using the specextract tool and created response 
files with the mkacisrmf and mkarf tools. We grouped each 
spectrum to have at least 50 counts per spectral channel. 



3 X-RAY SPECTRAL ANALYSIS 

We fit the data with XSPEC version 12.5. In, using the Sur- 
face Thermal Emission and Magnetospheric Scattering and 
assumed a reference gravitational redshift of 0.306 for the 
neutron star. We calculated the unabsorbed fluxes in the 0.5 
to 6.5 keV energy range using the cflux model in XSPEC. 
The distance of the SGR is assumed to be 48.1 kpc (Macri 
et al. 2006). Uncertainties reported throughout the paper 
correspond to 68% confidence limits of parameters, unless 
indicated otherwise. 

Uncertainties in the amount of interstellar absorption of 
X-ray photons often hampers efforts to precisely determine 
the intrinsic spectral shape of these sources. Furthermore, as 
a free parameter, hydrogen column density becomes strongly 
correlated with the determined temperature of the X-ray 
source essentially affecting the resulting best fit value (see, 
e.g. Durant & van Kerkwijk 2006 and Giiver et al. 2012 for 
a detailed discussion). Therefore, an accurate estimation of 
the hydrogen column density is crucial to deduce the spectral 
characteristics of the source. The fact that SGR 0526—66 is 
embedded in the supernova remnant N49 presents an advan- 
tage in this respect because it is possible to obtain an inde- 
pendent measurement of the hydrogen column density using 
the SNR. For that purpose, we extracted a SNR spectrum 
from the longest exposure and modeled it. Similar to the ear- 
lier studies of N49 (Park et al. 2003; Bilikova 2007; Tiengo 
et al. 2009; Park et al. 2012), we fit the remnant with a two 
component model, consisting of two plane-parallel shocked 
plasma functions (Borkowski et al. 2001, vpshock in XSPEC) 
both absorbed by the interstellar medium, assuming a solar 
abundance as given by Anders & Grevesse (1989). We ob- 
tained an adequate fit with a %t = 1-44 for 99 degrees of 
freedom. We note that the relatively high x 2 is mostly due a 
small number of individual energy bins, which could be emis- 
sion lines that can not be resolved by the ACIS-S detector. 
Our best fit parameters are: A<h = (0.15± 0.03) x 10 22 cm -2 , 
fcTi = 0.54 ± 0.01 keV, n > 3.6 x 10 13 s cm -3 , 
kT 2 = 1.11 ± 0.4 keV, r 2 = (1.21 ± 1.1) x 10 10 s cm -3 , 
where fcTi, kT 2 , n and t 2 denotes the temperature of the 
plasma and the upper limit on the ionziation timescale for 
each component. The elemental abundances we found are as 
follows: Ne/Ne = 0.90 ± 0.10, Mg/Mg = 0.64 ± 0.10, 
Si/Si = 0.66 ± 0.08, S/S = 1.01 ± 0.25, and 
Fe/Fe = 0.38 ± 0.04, all the other abundances were set 
to solar values. These results, especially the hydrogen col- 
umn density, are in general good agreement with the results 
of Park et al. (2012). We use the hydrogen column density 
we inferred from this fit as a fixed parameter when analyzing 
the spectra of SGR 0526—66 in the remainder of the paper. 

We simultaneously fit all six spectra with the STEMS 
model. STEMS model assumes a fully ionized and strongly 
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Figure 1. True color image of the supernova remnant N49. Red, green and blue corresponds to the 0.2-2.0, 2.0-4.0 and 4.0-10.0 keV 
ranges, respectively. Source (solid line) and SNR (dashed lines) extraction regions are also shown in cyan. 



magnetized hydrogen atmosphere on the surface of the neu- 
tron star, which determines the general spectral character- 
istics of its X-ray emission (Ozel 2001, 2003). In the mag- 
netosphere, these surface photons are further scattered by 
mildly relativistic charges (Lyutikov & Gavriil 2006). The 
model has four parameters: the effective temperature, mag- 
netic field strength at the surface as well as the resonant 
scattering optical depth and the velocity of charged parti- 
cles in the neutron star magnetosphere (see e.g., Giiver et al. 
2007, 2008, 2011 for details). First, we allowed all STEMS 
parameters to vary in all observations. Such a fit resulted 
in a x 2 /dof of 1.074 for 515 degrees of freedom (dof). Re- 
sulting effective surface temperature (~ 0.35 keV), magnetic 
field strength (~ 3.61 x 10 14 G), the magnetospheric scatter- 
ing optical depth (~ 5.66) and the average particle velocity 
(~ 0.55) did not show statistically significant or systematic 
variations between different observations. However, due to 
rather low signal to noise ratio of individual X-ray spectra, 
the error in individual model parameters were large. In order 
to obtain more constrained results on the effective temper- 
ature (therefore, the apparent emitting radius) and follow 
its time evolution, we linked the model parameters for mag- 
netic field strength, magnetospheric scattering optical depth 
and the average particle velocity. This way, we obtained a 
X 2 /dof of 1.086/530. We note here that, despite adding more 
dof to the fit, the change in the fit statistics was negligible, 
further indicating that the other parameters do not show 
statistically significant variations. We obtained the follow- 
ing best fit parameter values when only the surface effective 
temperature and model normalization were allowed to vary 
between observations: the surface magnetic field strength 
B= (3.74to'i4) x 10 14 G, magnetospheric optical depth to 
resonant scattering r = 5.37^'^, the average velocity of 
particles /3 = 0.52 ± 0.03. We then used the cflux model 
to calculate the unabsorbed X-ray flux in the 0.5—6.5 keV 
range for each observation. We present these flux measure- 
ments together with the effective surface temperature values 
for each observation in Table [2] We also show in Figure [2] 



the 68% confidence contours for the surface effective tem- 
perature and flux values. Both Figure [2] and Table [2] show 
that although the X-ray flux of SGR 0526—66 decreases by a 
factor of 20% over the last 9 years the inferred temperature 
values are fairly constant during this period. 

Given the fact that the surface temperature remains 
constant, we repeated our simultaneous fit by linking the 
surface effective temperature among all observations as well, 
in order to further constrain it and allowed the normaliza- 
tion, hence the emitting area of the STEMS model, to vary 
between observations. This way we obtained a xVdcf of 
1.083/535. Figure |3 shows the data together with the best 
fitting model and fit residuals. As expected, the best fit val- 
ues are very similar to what we found earlier: We obtain that 
the magnetic field strength at the surface is B= (3.73lg'Jg) x 
10 14 G, the optical depth to resonant scattering at the neu- 
tron star magnetosphere is r = 5.47^q 49 and the average 
velocity of the magnetospheric particles is P = 0.52 ± 0.03 
and an effective temperature of 0.3551q'q24 keV. We note 
that the best fit value of the surface temperature does not 
significantly depend on the fixed hydrogen column density. 
Even when we allow the hydrogen column density to be a 
free parameter in the fit to the spectra of SGR 0526—66, we 
obtain a value of iV H = (0.144 ± 0.012) x 10 22 cm -2 , which 
is consistent with the value we obtained from analyzing the 
spectrum of the surrounding supernova remnant. The un- 
absorbed flux and the emitting radius values obtained from 
this fit is given in Table [2] In Figure |4j we present y^/Aoi 
contours over the surface effective temperature and magnetic 
field strength. 



4 TIMING ANALYSIS 

We performed timing analysis using all available CXO 
observations to determine the spin period evolution of 
SGR 0526—66 and to uncover variations in the pulsed frac- 
tion. The time resolution of CXO observations in various 
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Table 2. X-ray spectral fit results for SGR 0526-66. 



Date Unabsorbed Flux" kT c ^ a Unabsorbed Flux 6 Radius 6 

(10~ 12 erg s" 1 cm" 2 ) (keV) (lCT 12 erg s" 1 cm" 2 ) (km) 



2000-01-04 


1.33±0.02 


n 37+0.04 

U - CW -0.04 


1.31±0.02 


13.33±1.56 


2001-08-31 


1.28±0.02 


36+ 04 


1.26±0.01 


13.09±1.54 


2009-07-18 


1.06±0.02 


36+ 05 


1.04±0.02 


11.88±1.37 


2009-07-31 


1.09±0.03 


n 3^+0.05 
u -' 3o -0.04 


1.07±0.02 


12.05±1.39 


2009-09-16 


1.04±0.02 


n q7 +0.03 
u -°'-0.05 


1.00±0.02 


11.71±1.35 


2009-09-19 


1.05±0.02 


n 34+0-04 


1.03±0.02 


11.85±1.37 



a Calculated from each individual data set when all the model parameters were linked between the observations, except for the effective 
temperature and normalizations. 

b Calculated from each individual data set when all the model parameters were linked between the observations, except for 
normalizations. 
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Figure 2. 68% confidence contours of unabsorbed 0.5 - 6.5 keV flux and surface effective temperature values for individual data sets 
when all other parameters were linked between different observations. 



sub-array modes is approximately 0.4 or 0.8 s, which is not 
ideal but sufficient for timing purposes. Using the same se- 
lection regions described in Section [3j we extracted source 
events for each observation and applied barycentric correc- 
tion using the axbary tool. 

We also included the deep XMM-Newton observation 
performed in 2007 (ObsID 0505310101) in our timing inves- 
tigation. We calibrated the EPIC-pn data using SAS v. 10.0.0 
and the calibration files as of October 2010. We extracted 
source events from a circular region with 10" radius. Note 
that the point spread function of XMM-Newton is not accu- 
rate enough to completely resolve the pulsar from the super- 
nova remnant; therefore, some fraction of unpulsed emission 
is expected to originate from the remnant (see Tiengo et 
al. 2009 for the details of contribution from the supernova 
remnant). We used the barycen tool of SAS to convert each 
event arrival times to that of the Solar system barycenter. 

To search for the pulsed signal from SGR 0526—66, 
we employed a technique (Buccheri et al. 1983) with 
the number of harmonics set to m = 2. We performed the 
search in a period range between 8.0 and 8.1 s. We detect the 
pulsed signal with high significance in the first three CXO 



data sets as well as in the XMM-Newton data, while the de- 
tections in the CXO data sets with observation IDs 10806, 
10807 and 10808 were marginal. In Table [3] we present only 
the statistically significant measurements of the spin period 
of SGR 0526—66, together with the chance probability and 
the Z 2 , power of each measurement. To determine the rate of 
change of the spin period, we fit the measured periods with a 
first order polynomial. We obtain a good fit with a spin down 
rate of (4.019 ± 0.494) xlO" 11 s s~\ We present the evolu- 
tion of the spin period as well as the best fit spin down rate 
of SGR 0526—66 in Figure [5] Note that our result regarding 
the spin evolution is based on only four measurements ob- 
tained over nine years. A more precise measurement of the 
spin down rate require more frequent deep observations. 

Finally, we constructed the pulse profile and calculated 
the root-mean-square pulsed fraction for each observation 
with significant spin period detection. The rms pulsed frac- 
tion (PF) is calculated as, 



PF = (^ ( £ (Ri " Ravc)2 ~ AR ' 2) j /Rave ' (1) 
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Energy [keV] 

Figure 3. All six X-ray spectra of SGR 0526—66 obtained since 2000 (crosses), best fitting model curves (solid lines) and the fit residuals 
(lower panel). 
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Figure 4. Variations of x 2 /dof as a function of surface effective temperature (left panel) and magnetic field strength (right panel) as 
inferred from the X-ray spectral fit. Dashed and dashed-dotted lines indicate the 68% and 90% confidence levels for each parameter. 



where N is the number of pulse phase bins (N=16), 
Ri is the source count rate in each phase bin, ARi is the 
associated uncertainty in the count rate, and Rave is the av- 
erage count rate of the pulse profile. We present the rms 
pulsed fraction values in Table [3] We find that rms pulsed 
fraction of SGR 0526—66 is very low and remains con- 
stant around 4% among CXO observations, while pulsed 
fraction obtained from XMM-Newton observation is signifi- 
cantly lower (~1.5%). Note, however, the fact that the rms 
pulsed fraction value is normalized by the average count rate 
of the pulse profile and any blended emission from the su- 
pernova remnant would increase the average rate and reduce 



the pulsed fraction. It is likely that the drop of rms pulsed 
fraction only seen in XMM-Newton observation is not in- 
trinsic to the source. 



5 DISCUSSION AND CONCLUSIONS 

We performed a systematic analysis of the archival Chan- 
dra observations of SGR 0526—66. Fitting the X-ray spec- 
tra with a strongly magnetized atmosphere model allowed 
us determine the strength of the surface magnetic field 
as B = 3.73 x 10 14 G and the effective temperature as 
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Table 3. Results of the timing analysis 



Obs ID 


Obs Date 
(MJD) 


pa 

spin 

(s) 


Maximum Z| 
Power 


Probability' 


Pulsed Fraction" 


747 


52415.400 


8.044(2) 


19.796 


5.5xl0~ 4 


0.038(10) 


1957 


52030.200 


8.0469(8) 


24.518 


6.3xl0~ 5 


0.042(9) 


XMM 


55044.000 


8.0546(7) 


19.027 


7.8xl0~ 4 


0.015(4) 


10123 


55090.500 


8.056(5) 


15.130 


4.5xl0" 3 


0.044(14) 



a 68% confidence limits of the measurements are in the last digit and given in parentheses. 

^Probability of obtaining the peak power in the Z| power spectrum if there was no periodic signal present. 

01/04/01 03/12/27 06/09/22 09/06/18 
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Figure 5. Spin period history of SGR 0526—66. 

0.355 keV. We also obtained the optical depth to resonant 
scattering in the magnetosphere, r=5.47, and the average 
velocity of the magnetospheric particles, 0.52c. All these 
parameters remain constant over the course of nine years 
within uncertainties. 

What has been observed to change over this interval is 
the source flux: it decreased by about 20% from 2000 to 2009 
(see also Park et al. 2012). Such a flux decay can result from 
a w 5-6% decrease in the surface temperature. This varia- 
tion is comparable to the uncertainties in individual tem- 
perature measurements at the 1-a level. Therefore, we can 
not unambiguously rule out a variation in the temperature. 
Nevertheless, no systematic variation in the best fit values 
of the surface temperature has been observed over the nine 
years. Therefore, it is likely that the surface temperature of 
SGR 0526—66 remained constant throughout this period. 

The decline in the observed flux can also be achieved 
with a decrease in the radius of the emitting region by about 
10%. We obtain through our spectral fits that the radius was 
about 13.3 km in 2000/2001, and about 11.7 km in 2009, as- 
suming a distance 48.1 kpc (Macri et al. 2006). It is possible 
that the thermally emitting hot spot on the surface of the 
neutron star might still carry the effects of burst induced 
heating in the late 70s and early 80s and dissipates over 
time without going through a major change in its temper- 
ature. It is important to note that the large radius values 
we deduce from the spectral fits are in agreement with the 
very low rms pulsed fractions we inferred from the timing 
analysis. 




We determined the spin period history of SGR 0526—66, 
which is in accord with the results of Tiengo et al. (2009) and 
derived the period derivative as 4.02 x 10~ 11 s s~ 1 . Assuming 
that the spindown of the neutron star is due to a magnetic 
dipole braking, the inferred dipole magnetic field strength of 
the neutron star is calculated as ~ 5.7x 10 14 G, which agrees 
with the surface magnetic field strength as found from the 
X-ray spectra using the STEMS model. Making the same 
spindown assumption, we calculate the characteristic age of 
the pulsar as r s( j = P/2P ~ 3200 yr. This is comparable 
to, but somewhat lower than, the estimated Sedov age of 
the supernova remnant, which is 4800 yr (Park et al. 2012). 
Because SGRs and AXPs exhibit episodic variable spindown 
behavior (see, e.g., Woods et al. 2002), the characteristic age 
should only be taken as a rough indicator of the true age of 
SGR 0526-66. If SGR 0526-66 is indeed associated with 
N49, it requires a rather high radial velocity for the neutron 
star compared to some other AXPs and SGRs (Gaensler et 
al. 2001; Kaplan et al. 2009). 

Long term observations of magnetars may play an im- 
portant role in our understanding of the magnetic field decay 
and its effects on the neutron star crust and long term cool- 
ing. However, observations of magnetars over the last cou- 
ple of decades showed that these neutron stars are rarely 
in quiescent state: They show giant flares (Hurley et al. 
1999; 2005), outbursts (Ng et al. 2011), glitches (Dib et 
al. 2007), which likely affect the long-term cooling prevent 
any observational constrains on the theoretical calculations. 
Although located at a large distance and residing within a 
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SNR, SGR 0526—66 is an exception in this regard. Since the 
giant flare observed in 1979 and the bursting activity that 
lasted till 1983, no burst activity has been reported from 
this source, enabling us to assume that the source has been 
steadily cooling in quiescence. This is also confirmed by our 
spectral results presented in Section [3] One other advantage 
of this source is its possible association with the supernova 
remnant N49, providing a reasonable age estimate for the 
neutron star. 

Aguilera et al. (2008a) calculated cooling curves for 
strongly magnetized neutron stars by taking into account 
the Ohmic dissipation (rohm) and Hall drift (rnaii) pro- 
cesses, where the latter is employed for a rapid decay during 
early stages (see eq. 17 in Aguilera et al. 2008a). We can 
compare our results with the magnetic neutron star cooling 
curves calculated by Aguilera et al. (2008a, b) by using a 
conservative age range for SGR 0526-66, together with the 
surface temperature that we determined from six Chandra 
observations covering nearly a 10 year time interval. We take 
the age of the neutron star to be 3000 — 7000 yr, which spans 
the spindown age of SGR 0526-66 and the supernova rem- 
nant age reported by Park et al. (2003, 2012). In Figure [(J 
we present the cooling curves for magnetars with various 
initial magnetic field strengths, along with the current sur- 
face temperature of SGR 0526—66. Our comparison shows 
the cooling trend of SGR 0526—66 matches the theoretical 
cooling curves corresponding to two sets of parameters, an 
initial magnetic field strength of 10 16 G, Tohm of 10 6 yr and 
TOhm/THaii = 5 x 10 2 yr, or the curve that corresponds to an 
initial magnetic field strength of 1-5 x 10 15 G, Tohm of 10 6 yr 
and TOhm/THaii = 5xl0 3 yr. 

On either cooling track, SGR 0526—66 is expected to 
enter a phase of more rapid cooling in the near future, based 
on its age and its current temperature, which we reported 
here. If its temperature indeed begins to decay rapidly, ob- 
servations of the source over the next decade with planned 
X-ray telescopes, such as Advanced Telescope for High En- 
ergy Astrophysics, may reveal a more significant drop in the 
spectral temperature, accompanied by a further flux decay, 
as long as the source remains in its deep quiescence. 
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